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• We report on the status of our simulations of two-colour QCD with two flavours of Wilson 

quarks, at nonzero chemical potential and diquark source. Preliminary results are presented 
' for the static quark potential, gluon propagator, and meson and diquark spectrum. 
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O ■ §1. Motivation 

o ■ 
<N : 

T— I , QCD at non-zero density is currently an area of great activity, as the wide range 

' of presentations at this workshop bears witness. What is missing at present, how- 

. ever, is a first-principles, nonperturbative approach to QCD at high density and low 

, temperature, without which we have little hope of obtaining quantitative predictions 

I ' regarding the expected complex phase structure in this regime. Lattice QCD would 

be the natural candidate for such an approach, but it is severely hampered by the 
r~| ! fact that the euclidean action becomes complex when a non-zero chemical potential 

■ is introduced. Ref. gives an overview of lattice QCD at non-zero density and the 
• '-j . problems encountered. 

I We may avoid these problems by studying model theories which share some of 

' the salient features of QCD, while having a positive definite fermion determinant also 

for non-zero chemical potential. Examples of such theories are 2-colour QCD,'^ QCD 
with adjoint fermions,"^' NJL models,^ and QCD with isospin chemical potential.!^ 
Here we will be concentrating on 2-colour QCD, which in the gauge sector is expected 
to be very similar to real QCD, although the baryon sector and chiral properties of 
the two theories are radically different. 

It is of fundamental importance that a theory to be studied using lattice regu- 
larisation have a well-defined continuum limit. In the context of non-abelian gauge 
theories, this implies that the theory should exhibit asymptotic freedom. For SU(2) 
the second term in the perturbative /^-function changes sign for Nf > 5.6, with rad- 
ical consequences for the infrared behaviour of the theory; nonperturbative effects 
reducing this critical value of A'^^ still further cannot be ruled out. Consequently 
we prefer to work with Nf = 2 so that both asymptotic freedom (and hence a 
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continuum limit) and confining behavour in the infrared are guaranteed. Most sim- 
ulations up to now have been performed using staggered fermions, where either A'^^ 
is a multiple of 4, or a square-root of the fermion determinant is taken, making the 
action non-local. To avoid the potential pitfalls this creates we are here using the 
Wilson action with Nj = 2 quark flavours. This action has also been employed by 
the Hiroshima groupP'''^ albeit with a different gauge action, so our results may be 
directly compared to theirs. 

It is also worth pointing out that for k = Kc Wilson fermions, unlike staggered 
fermions, have the same global symmetry properties as continuum 2-color QCD, 
as encoded by the Dyson index (3 = 1^ which has consequences for the physical 
spectrum of the model as analysed in Ref . !^ . 

We now list some physical motivations for studying this model. 

1.1. Phase diagram 

There are indications, from studies at nonzero isospin-chemical potential, that 
the curvature of the phase transition line in the (T, /i)-plane for low /i is the same as 
in real QCD. If this can be confirmed, we may gain insight into the phase diagram 
of real QCD by studying 2-colour QCD, where numerical studies are not restricted 
to small /u/T. 

At the other end of the phase diagram, at low T and high /i, a possible tricritical 
point has been identified.^ A confirmation of this would be welcome. Another 
interesting issue at low T and intermediate /i is whether 2-colour QCD has a normal 
matter phase, i.e. a phase with nonzero baryon number density and zero diquark 
condensate — or whether the finite-// transition is directly from the vacuum to 
a superfluid phase. The available evidence, as well as chiral perturbation theory 
studies, indicate the latter, but the matter has not yet been finally settled. 

1.2. Spectrum 

Previous studies of 2-colour QCD with Wilson fermion^''^ have indicated that 
the vector meson becomes lighter with increasing chemical potential, in line with 
the predictions of Brown-Rho scaling^^' which have been put forward as a possible 
explanation for the CERES results.^ Such an in-medium effect, if confirmed, would 
provide a new testing ground in heavy-ion experiments such as the one planned at 
GSI. Closely related to in-medium modification of vector meson properties is the 
possibility of vector condensation, which may occur instead of or in addition to 
scalar diquark condensation at high density.^'^ 

The behaviour of the pseudoscalar isosinglet diquark at high density is also of 
interest, since it shares the quantum numbers of the r]' meson but is accessible via 
standard lattice spectroscopic techniques. Its becoming light as fi increases may 
signal the restoration of U(l)^ symmetry in dense matter.C^ 

Finally, we hope to compare the actual behaviour of the pion and the C-even 
and -odd scalar diquark masses, as well as the baryon density and chiral and diquark 
condensates, with the predictions from chiral perturbation theory .'Si The agreement 
(for the chiral condensate in particular) is striking in the staggered adjoint modeP^ 
which is in the same global symmetry class (/? = 1) as the model we are studying 
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here. At some point, of course, chiral perturbation should cease to be accurate once 
/i becomes comparable with the mass of the vector diquark, and it is important to 
know the range of parameters where this breakdown can be anticipated. 

1.3. Gluodynamics 

Gluodynamics is arguably the sector where results from two-colour QCD have 
most direct relevance to real QCD. At zero chemical potential, the gluodynamics 
of SU(2) and SU(3) are so similar that SU(2) simulations have often been used 
as a 'cheaper' way of obtaining results relevant to QCD. This is also where we 
most directly can study effects of fermion loops, including Pauli blocking at high 
density. In contrast, all numerical results for spectrum and condensates to date 
might conceivably have been obtained in the quenched approximation. 

Among the quantities of particular interest are: 

• Polyakov line screening and static quark potential. Previous simulations of 
SU(2) with Wilson fermions have reported signs of a deconfinement transition 
at large while no such transition has been observed with staggered fermions. 

• Gluon propagator. This quantity is an essential ingredient also in the Dyson- 
Schwinger equation approach to gauge theories .^'^ These approaches can be 
applied equally to SU(2) and SU(3) gauge theories, so input from SU(2) lattice 
simulations may provide a valuable check on the methods and assumptions 
used. 

• Other quantities that may provide further insight into the phase transition 
include glueball massefl^ and the gluon condensate.!^ 



§2. Algorithm 



Above the critical chemical potential Hc for diquark condensation, it is essential 
to have a diquark source. This serves two purposes: 

1. It lifts (near-)zero eigenvalues from the Goldstone modes, improving the condi- 
tion number and avoiding critical slowing down. 

2. It also enables the determination of diquark condensates and other anomalous 
propagation. 

The action, including diquark source terms, can be written 



+ J 



(2-1) 



where j,j are sources for the diquark and antidiquark respectively, and M(^) is the 
usual Wilson-fermion matrix. If we introduce M{ijl) = ^^M{n) and the new field 
variables 



Xi = ^1 , Xi = V'Ts ■ 
the action can be rewritten as 



92 



T2CiIjX , 



s 




Xi 
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XQ[f,,j,j]X. 



(2-2) 



(2-3) 
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The matrix Q is hermitean, and the determinant is negative definite if the two 
sources are equal: j = j, so this representation lends itself to simulation using 
standard algorithms. 

We will use the two-step multibosonic algorithm!^ which has been shown to be 
more efficient than Hybrid Monte Carlo (HMC) in the dense phasej^ and which is 
ergodic even in the presence of singularities. However, in the initial phase, we have 
also performed simulations using HMC without a diquark source, and these are the 
results we will quote in the following. 

§3. Results 

In order to establish a baseline, we have first performed simulations on 8'^ x 16 
lattices at = 0, for several values of /? and k, using the standard Wilson action 
both for gauge fields and for fermions. The parameters are given in table along 
with the associated average spatial and timelike plaquette values. We have sampled 
configurations every 4 trajectories. 
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dt 




acc 




□t 


Name 


1.7 


0.1780 


0.0125 


7152 


84% 


0.47386(6) 


0.47367(6) 


coarse (c) 


1.8 


0.1725 


0.02 


3856 


70% 


0.50007(8) 


0.49983(8) 






0.1740 


0.01 


216 


89% 












0.0125 


272 


80% 


0.5056(2) 


0.5052(2) 






0.1750 


0.01 


600 


80% 


0.5140(3) 


0.5139(3) 


light (1) 


1.9 


0.1650 


0.02 


2316 


83% 


0.52124(10) 


0.52098(10) 


fine (f) 



Table I. Hybrid Monte Carlo simulation parameters, fi — 0. dt is the HMC timestep, while 'acc' 
denotes the acceptance rate. Ds and are the average spatial and timelike plaquette values 
respectively. 

3.1. Setting the scale 

We have attempted to set the scale from the string tension by computing Creutz 
ratios. Figure ^ shows the results for the coarse lattice (/3 = 1.8, ac = 0.178) which 
is where we have the best statistics. It is evident from the figure that our statistics 
are still not sufficient to determine the string tension from Creutz ratios with any 
degree of confidence, and from the other parameter values the situation is even worse. 
Based on the trend displayed by the data, we tentatively estimate a^a ~ 0.37(5). 
Taking the string tension to be ^/a ~ 420 MeV, this gives a lattice spacing a ~ 0.29 
fm or « 690 MeV. 

In Fig. 121 we show Creutz ratios for R = 2 only for all our parameter values. 
We also include the i? = 3,4 data for the light lattice (/3 = 1.8, k = 0.175), to give 
an indication of the trend. We may attempt to determine ratios of lattice spacing 
by simply taking the ratios of Creutz ratios at {R = 2,T = 3) where the statistics 
permit this. The results are shown in table ^ These numbers are likely to be 
overestimates, since ratios for R = 2 appear to increase rather than decrease with 
T. 

We have also attempted to determine the lattice spacing ratios by matching 
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Fig. 1. Creutz ratios for the coarse lattice (/3 = 1.7, k = 0.178). 
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Fig. 2. Creutz ratios for R = 2, all k) values. 







Creutz 




Gluon 






K 


Ra 


Ra 


Rz 




1.7 


0.1780 


1 


1 


1 




1.8 


0.1725 


0.95 


0.84(1) 


1.06 


9.7 




0.1740 


0.93 


0.76(2) 


1.06 


2.0 




0.1750 


0.89 


0.69(3) 


1.06 


2.7 


1.9 


0.1650 


0.92 


0.74(1) 


1.14 


20 



Table II. Ratios of lattice spacings Ra = a/flc, where is the lattice spacing on the coarse lat- 
tice. 'Creutz' denotes the ratios of Creutz ratios for J? = 2, T = 3, while 'Gluon' denotes the 
outcome of matching the gluon propagator as described in Ref. I22L with Rz being the ratio of 
renormalisation constants Za/Z^. 
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the (Landau-gauge) gluon propagator as described in Ref. [22]) . We have imposed 
a cyhnder cut, selecting points within a radius 1.2 x 7r/8a from the diagonal in 
momentum space, and varied the ratios of lattice spacing a and gluon renormalisation 
constant Z3 to match the gluon propagator for the different parameters. In Fig. |21 
we show the gluon form factor before and after matching. The resulting a and Z3 
ratios are given in table Ull 



A |i=1.7 K=0.1780 - 

o p=1.8 ic=0.1725 - 

□ ic=0.1740 ~ 

, ic=0.1750 - 

P=L9 K=0.1650 I 



A P=l 


7 K=0.1780 


p=l 


8 K=0.1725 


□ P^l 


8 K=0.1740 




8 K=0.1750 




9 K=0.1650 



Fig. 3. The gluon propagator form factor in Landau gauge, as a function of momentum. Left: 
before matching of lattice spacings; right: after matching. 



This procedure should be treated with some caution, since it can fail on three 
counts. Firstly, in the presence of dynamical fermions the renormalised gluon prop- 
agator will depend on the renormalised fermion masses, which vary between our pa- 
rameter sets. This may explain the discrepancy between the data for {(5 = 1.8, k = 
0.175) and (/3 = 1.9, k = 0.165) at large momenta since, as we shall see, these 
two have the lightest and the heaviest quark mass respectively. Secondly, on these 
coarse lattices it is likely that we are not in the scaling regime so the renormalised 
propagator will have a substantial cutoff dependence. Thirdly, since the volume is 
fixed in lattice units, these lattices will have different physical volumes and finite 
volume effects should cause a slight mismatch in the infrared, where our procedure 
has left them almost perfectly matched. In particular, for the lightest quark mass, 
corresponding to the finest lattice, we see quite large lattice artefacts at intermediate 
momenta, which may arise from finite volume effects. For the other lattices however, 
the physical volumes are so large that such effects are likely to be quite small. 

With these caveats in mind, the two methods may complement each other and 
give an idea of the systematic uncertainties in the lattice spacings. It would clearly 
be preferable to compute the static quark potential using smeared Wilson loops, but 
even this may fail to provide a reliable lattice spacing on these small lattices. 

3.2. Meson and diquark spectrum 

We have computed correlators, using point sources, for the pseudoscalar (vr), 
scalar {5) and vector (p) mesons, as well as for the scalar, pseudoscalar and vector 
diquarks. For the vector diquark we have not been able to extract a mass, while the 
scalar and pseudoscalar diquarks are exactly degenerate with the vr and 6, as they 
should be at ^ = 0. In the following, we will therefore only quote our results for the 
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meson spectrum, with the main quantity of interest being mT^/mp. 
In Figure 0] we show the effective 



meson masses for (/? = 1.7, k = 0.178). 
We see, firstly, that with point sources 
on these lattices, we do not get a good 
plateau for either of these correlators, 
and secondly, that the 5 correlator is 
still too noisy to determine the mass 
with any confidence, although clearly 
ms ^ rriTr, implying we are in a phase 
where chiral symmetry is spontaneously 
broken. Although we do not see a 




plateau, we can still fit the vr and p cor- 
relators to a double-exponential (cosh) Fig- 4. Effective masses for (/3 = 0.17, k = 
with a good to reasonable . The re- 0.178). 
suiting masses are given in table II I II 

However, as the absence of a plateau in the effective-mass plot would indicate, the 
fit values are not stable, and one should add a systematic uncertainty of 1-2% in all 
quantities from the variation with the fit range. 





ft 


m,r 




rup 




range 


m-^/nip 


1.7 


0.1780 


0.800(2) 


0.97 


0.870(3) 


1.4 


3-5 


0.920(3) 


1.8 


0.1725 


0.807(3) 


1.3 


0.859(4) 


0.82 


4-6 


0.939(4) 




0.1740 


0.731(9) 


2.5 


0.820(14) 


1.6 


3-5 


0.892(15) 




0.1750 


0.642(9) 


0.64 


0.726(15) 


1.4 


3-5 


0.885(16) 


1.9 


0.1650 


0.925(4) 


0.27 


0.973(5) 


0.12 


4-6 


0.951(4) 



Table III. Meson masses, m-n- and nip axe obtained by fitting the correlators to a cosh function 
on the timeslices denoted by 'range' (and their mirror images around the centre of the lattice). 
The errors are purely statistical bootstrap errors. 



Comparing tables Ulland lllll it appears that the lattice spacing changes somewhat 
more rapidly with k at fixed (3 than the vr-to-/? mass ratio. A more careful analysis is 
needed to identify the true lattice spacing and spectrum before it can be determined 
whether this is a real effect or merely a sign of shortcomings in our analysis. 

Finally, in Fig. [5] we show the scalar meson and vector diquark correlators as a 
function of time. Both correlators change sign, so we plot the absolute value in order 
to show them on a logarithmic scale. 



§4. Outlook 



We have performed initial simulations to explore the parameter space for 2- 
colour QCD with Wilson fermions. The next step will be to extend the simulations 
at one or more of these points to non-zero chemical potential. Preliminary studies 
show that HMC experiences a dramatic slowing-down as we approach the expected 
transition point at /U = m-,^/2. We will therefore be using the TSMB algorithm to 
explore densities beyond this point. 
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Fig. 5. Absolute value of the scalar meson (left) and vector diquark (right) correlators. 
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